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Abstract 

We study the dynamics of a quantum particle moving in a plane 
under the influence of a constant magnetic field and driven by a slowly 
time-dependent singular flux tube through a puncture. The known 
adiabatic results do not cover these models as the Hamiltonian has 
time dependent domain. We give a meaning to the propagator and 
prove an adiabatic theorem. To this end we introduce and develop 
the new notion of a propagator weakly associated to a time-dependent 
Hamiltonian. 



1 Introduction 

The model under consideration originates from Laughlin's jlSj and Halperin's 
[HI discussion of the Integer Quantum Hall effect. In the mathematical physics 
literature Bellissard [H] and Avron, Seiler, Simon [3] used an adiabatic limit 
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of the model (with additional randomness) to introduce indices. The indices 
explain the quantization of charge transport observed in the experiments |12j . 

In this paper we discuss some mathematical aspects of the existence of 
the propagator and the validity of the adiabatic approximation and propose 
how to overcome the difficulties originating from the strong singularity of the 
external field. 

Let us specify the model, summarize our results and introduce the nota- 
tion. The configuration space is M 2 \ {(0, 0)} and the model is considered in 
polar coordinates (r, 6). The vector potential A is the sum of a part for the 
homogeneous magnetic field of strength B > 0, 

B Br 2 
— {x x dx 2 - x 2 dxi) = — - d8, 

plus a part describing the flux <3> which varies in time, 

$ 1 $ 
— -p^j (xidx 2 - x 2 dxi) = — d9; 

Z7T \X\ Z 27T 

the real- valued function $ is assumed to be monotonous and C 2 . With the 
metric coefficients gn = 1, g 22 = r 2 , f? 12 = 0, the differential expression of 
the Hamiltonian acting in L 2 (IR + x [0,2ii[,rdrd9) is 

_ _ Jgc? k (-iW k - -A k ) 

= HL (A d rd +- (-id ----- ®\\ 
2m I r r r r 2 \ 8 he 2 he J J 

Our purpose is to study the response of the system if flux quanta hc/e are 
added adiabatically, i.e. the flux function is of the form t \— > $(t/r) with the 
time t varying in [0, r] for some r 3> 1. 

In a first step we analyze the case when $ is linear. Furthermore, we fix 
an angular momentum sector defined by —idoe ime = me im6 {m e Z), and use 
a slow time s, i.e.: the substitution s = —m + e/(hc)^(t/r). Also we are not 
interested here in keeping track of the behavior in the physical parameters e, 
h, c, 2m, so we set them all equal to one. This is our motivation to consider 
the operator 

1 1 / Br 2 \ 2 

H(s) = — d r rd r + - [s + — - in L 2 (R + , rdr). (1) 
r r z \ 2 ) 

In a second step we shall then show that our analysis generalizes to Hamil- 
tonians of the form H^(s)^j where C G C 2 is a monotone function. 
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H(s) is essentially selfadjoint on Cg°(]0, oof) iff s 2 > 1 []3|. For < s 2 < 1 
we impose the regular boundary condition as r — * 0+ (i.e.: a wavefunction 
belongs to the domain if it has no part proportional to the (square integrable) 
singularity r~' s '). This is in fact the most common choice, see [H] for a 
detailed discussion. The case s = is particular since the singularity in 
question is logarithmic but otherwise the situation is similar, see [T]. The 
Hamiltonian H(s) is unambiguously determined by specifying a complete set 
of eigenfunctions with corresponding eigenvalues, see below. 

The dynamics of the model should be defined by 

id s U T (s, s )ip = tH(s)U t (s, s )ip, U T (s , s )ip = ip, (2) 

where U T is unitary and ip is an arbitrary initial condition from the domain 
of H(sq). The existence of a propagator in this sense is, however, uncertain. 
The problem arises from the fact that the domain of H(s) is not constant in s, 
respectively that H(s) is not relatively bounded with respect to H(s). Thus 
the usual theorems which assure the existence of the propagator [14] and the 
validity of the adiabatic approximation [H E] are not directly applicable. 

A convenient way to see this is to consider the eigenfunctions. The oper- 
ator H(s) has a simple discrete spectrum; the eigenvalues are 

A n (s) =B(s + \s\ +2n + l), n6 {0,1,2,...}, (3) 

with the corresponding normalized eigenfunctions 

^( S ;r) = Cn ( S )rHL(N)(^)exp(-^ 

where 

/ b n(M+d/j/ 2 „| y/2 

c " (s) = U ) U» + w + dJ 

(\s\) 

are the normalization constants and L„ are the generalized Laguerre poly- 
nomials (see, for example, \S\). 
The derivative of H(s) equals 

ff( S ) = | + B. 

Notice that if \s\ < 1 then (p n (s) cannot belong to the domain Dom if(s) 
since H(s)ip n (s) ~ r~ 2+ ^ for r — > 0+. This means that H(s) is not relatively 
bounded with respect to H(s). 
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Remark that, on the other hand, the quadratic expression 

/■oo 

/ ip m (s;r) H(s)ip n (s;r)rdr 
Jo 

makes good sense. In order to avoid a complicated notation we shall denote 
it by the symbol ((p m (s), H(s)<p n (s)) even though the symbol cannot be taken 
literally and is therefore somewhat misleading. Furthermore, the derivative 
of the eigenfunction, <f n (s), belongs to L 2 (M + , r dr). Since the eigenfunctions 
are chosen to be real-valued it holds true that 

(ip n (s),ip n (s)) = 0. 

Let us also note that, similarly, if \s\ < 1 and s 2 ^ s' 2 then the eigen- 
function ip n (s) cannot belong to Dom H(s'). It is so because (as formal 
expressions) H(s') — H(s) = (s' 2 — s 2 )/r 2 + B(s f — s) and H(s')ip n (s;r) has 
a non-integrable singularity at r = 0. Hence Dom H (s) depends on s. 

It turns out that, following the strategy of Born and Fock [2j, the problems 
of existence and adiabatic approximation can both be handled: 

denote the eigenprojector onto C<p n (s) by -P n (s); it is differentiate as a 
bounded operator. The hard part of our work consists in showing that 

00 

is a bounded operator. This is stated in Lemma El It requires work because 
its matrix elements have bad off-diagonal decay, see Lemma |H (which is 
formulated for the unitarily equivalent operator Q). 
Now 

00 

H AD {s) := H(s) + - V Pn(s)P n (s) 

T * — ' 

n=0 

has a propagator which is well defined in the usual way, i.e. 

id s U AD (s, s )ip = tH ad (s)U ad (s, s )ip, U AD (s , s )ip = ip, (4) 

for 1/} G Dom (H AD (s )). To see this notice that U A d can be computed by its 
action on the eigenbasis: 

tt I \ I \ —i T f S A n (u)d« / \ 

U AD {s,S )ip n {s ) = 6 J °a ¥n{s)- 

Furthermore, X n (s) — A n (0) is bounded in n and so U A d(s, s ) DomH AD (s ) = 
Dom H AD (s). 
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Since H(s) — H AD (s) is bounded the domains of H(s) and H AD (s) are 
identical. By time-dependent transformation a natural candidate for the 
propagator of H(s) is 

U T (s, s ) := U AD (s, 0)C(s, s o )U AD {0, s ) (5) 

where C(s, sq) is defined by 

id s C(s, s ) = -Q T (s)C(s, s ), C(s , s ) = I, (6) 

with 

Q T (s) :=U AD (0,s) (ijrP k (s)P k (s)) U AD (s,0). (7) 
V k=o J 

Since ||Q r (s)|| is locally bounded the propagator C(s,s ) is well defined by 
the Dyson formula. 

The adiabatic approximation problem is settled in Proposition fTTI were it 
is shown that 

\\U T (s,0)-U AD (s,0)\\=O 

It remains unclear, however, whether C(s,sq) preserves the domain of 
H(0) and therefore whether the propagator U T (s,s ) is actually related to 
the Hamiltonian H(s) in the usual sense. To handle this problem we develop 
the general concept of weak association of a propagator and a time dependent 
Hamiltonian. We can show that U T is weakly associated to H(s) and that 
the Schrodinger equation (J2J) is fulfilled in the sense of distributions. 

We shall use the following notation. The symbol V(s) stands for the 
unitary operator which sends all eigenstates at time to the corresponding 
eigenstates at time s, i.e. 

V(s)(p n (0) = cp n (s) Vn G Z + (8) 

(here and everywhere in what follows Z + stands for the set of nonnegative 
integers). Further set 

oo 

W(s) = KOO^if (s)V(s) = A «( s ) ( 9 ) 

n=0 

and 

oo 

Q(s) = ^ n (s)P n (0) (10) 

n=0 
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where 



^n(s) = / X n (u) du. 

Jo 

Remark that the adiabatic propagator decomposes as 



U AD (s,s ) = V(s)e- i ^- u ^V(s )- 



The paper is organized as follows. In Sections [2] and El we do the analysis 
necessary to prove the boundedness result stated in Lemma El Section 0] is 
devoted to the existence problem for the propagator. In Section El we prove 
the adiabatic theorem in Proposition [TTJ The result is then extended to a 
more general time-dependence in Sectional 

A rather independent part of the paper is the Appendix where we propose 
the notion of a propagator weakly associated to a time-dependent Hamilto- 
nian. We indicate cases where the weak association can be verified while 
the usual relationship between a propagator and a Hamiltonian is unclear or 
even is not valid. In particular, this concept was inspired by the situation we 
encountered in the present model. We believe, however, that this idea need 
not be restricted to this case only and that it might turn out to be useful in 
resolving this type of difficulties in other models as well. 

2 Auxiliary estimates of matrix operators 

Here we derive some auxiliary estimates that will be useful later when veri- 
fying assumptions of the adiabatic theorem. 

Lemma 1. Let A(o~) be an operator in l 2 (N) depending on a parameter o > 
whose matrix entries in the standard basis equal 



Then A(o~) is bounded, uniformly in a, and its norm satisfies the estimate 



Proof. The proof will be done in several steps. 

(i) Let K(a) be an integral operator acting in L 2 (M + , dx) with the integral 
kernel 




for m = n 
for m < n 
for m > n 



\\A{a)\\ < 24. 
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Let us show that 

\\K(a)\\ = —?—. 
11 V JU 2a + 1 

First we apply the unitary transform 

U : L 2 (R + , dx) -> L 2 (R, dy), U^(y) = e y/2 ^(e y ). (11) 
The inverse transform reads U~ 1, \p{x) = x -1 / 2 ^(rna;). Set 

K(a) = UK(a)U-\ 
One finds that K(a) is again an integral operator with the integral kernel 
K a {y,z)=isgn(y-z) e -("+V2)|»-*|. 

Hence K{a) is a convolution operator and it is therefore diagonalizable with 
the aid of the Fourier transform T on R. This means that 

where 



q( z ) = f e lzy sgn(y) e^ a+1/2)lyl dy 



2iz 



2 ■ 2 
Z 2 



It follows that 

\\K(a)\\ = \\FK(a)F- 1 \\ = \\q\\ 00 = ^- T . (12) 

a + 2 

(ii) Suppose that {ip}^ =l is an orthogonal system in L 2 (R + , dx) such that 

Vm,n G N, (il> m ,K(a)il> n ) = A(a) mn 

and 

\/n G N, H^nf = « > 0. 

Let P + be the orthogonal projector onto spanj^}^ in L 2 {R + ,dx). Then 
one can identify P + K(a)P + with K^ 1 A(a). Hence 

\\A(a)\\ = K\\P + K(a)P + \\<K\\K(a)\\. (13) 

fmj We shall construct an orthogonal system {il) n }™=\ described in the 
preceding point as follows. Consider the natural embedding L 2 ([n, n+l], dx) C 
L 2 (R+, dx), n G N. We seek *0n G L 2 ([n, n + 1], dx) in the form 

1p n = a n U n + (3 n V n + f n 
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where a n , (3 n G R, f n G L 2 ([n,n + l],ds), 

u n {x) = x a , v n (x) = x _fT_1 for x G [n, n + 1], 
and f n -Lu ni f n -Lv n . Suppose for definiteness that m < n. Then 

rm.+ l rn+1 

(ip m ,K(a)i() n ) = dx dyK, a (x 7 y)ip m (x)ip n (y) 



m J n 



= -i (u m , ip m ) (v n , ij) n ) . 

Furthermore, 

(tfj n ,K(a)tp n ) = / K c {x,y)ip n {x) ip n {y) dxdy = 

J n J n 

since K a (x,y) is antisymmetric, JC a (y,x) = —fC a (x,y). Consequently, it 
suffices to choose the real coefficients a n , (3 n so that 



Vn G N, (u n ,ip n ) = n a , (v n ,^ n ) = n 



-CT-l 



This system has a unique solution (a n , f3 n ). The function f n can be arbitrary. 
Its only purpose is to adjust the norms of the functions ip n so that they are 
all equal. Set 

N n {a) = \\a n u n + f3 n v n \\ 2 = / (a n x a + (3 n x~ a ~ l ) dx 

J n 

and 

k(o) = supA^(cr). 

neN 

One can choose the orthogonal system {i/Jn}n°=i so that H^/Vill 2 = k{o) for all 
n. According to (|T2|) and (fTT?|) we have 

IWII<|^. ("I 
(iv) It remains to find an upper bound on k(<j). Set 

£ n = n a , r] n = n~ a ~ 1 . 
Simple algebraic manipulations yield 

) ZnVn + (u n ,U n ) Vn 

™n(. cr ) = 2 ' 

) - (U n , V n ) 



Here 



(u n ,v n ) = ln(l + - 

V n , 

(u n ,u n ) = -i— ((n + l) 2ff+1 - ^ 2CT+1 ) , 

(Vn> = ^ViK 2<T_1 -^+ 1 )^ 1 )- 



w=(o-+-)ln(l + -). 



Set 



One can rewrite the expression for N n (a) as follows, 

2cr + 1 sinh(w) cosh(w) — w 

N n (a) = . 2 . 

n sinn (w) — ur 

Using an elementary analysis one can show that 

sinh(w) coshfu;) — w smh(w) cosh(w) — w , , . 

9/ x < — —7— T-7 — r < 4 cotgh (w) 

sinh (it?) — sinh(w) (sinh(w) — w) 

Hence 



NJo) < — ~ - < 12 ( 2(J + !)• 



4(2* +1) (l + tf^ + l 
Consequently, 

k(o) < 12(2(7+1). (15) 
From JUJ and jEJ it follows that ||A(a)|| < 24. □ 

Lemma 2. Let A(er) be an operator in l 2 (N) whose matrix entries in the 
standard basis equal 

for m = n 

A{a) mn ={ -if*(%) form<n 



■/*(£) form> 



n 



where 

1 — u a 

fa{u) = Y^"' U G ]°' 1 [' 

and a G [0, 1] is a parameter. Then A(o~) is bounded and its norm satisfies 
the estimate 

ll^)l!<| J + 4|A. 
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Proof. The proof will be done in several steps. 

(i) Let K(p) be an integral operator acting in L 2 (R + , dx) with the integral 
kernel 

, for*<</ 

/C a (x,y) 

Let us show that 



i for:r>y 



11^)11 <ttV. (16) 

This step is quite analogous to the proof of point (i) in Lemma HJ First we 
apply the unitary transform U defined in (|11|> . Set 



Kip) = UK{a)U-\ 
One finds that K(a) is again an integral operator with the integral kernel 

K a {y,z) =isgn{y - z) f^e^) e~^ 2 . 

Thus K{p) is a convolution operator which is diagonalizable with the aid 
of the Fourier transform T on R. This means that (j r K(a)J-'~ 1 ip) (z) = 
q(z)i/j(z) where 

q(z) = [ e^sgn(y)f a (e-\y\)e-\y\/ 2 dy. 
Jr 

A standard estimate yields 
It follows that 

\\K(o-)\\ = \\FK{o-)F- l \\ = \\q\\ 00 <iT 2 o-. 

(ii) Let Xn{ x ) De the characteristic function of the interval }n,n + 1[. 
The linear mapping 

oo 

J : / 2 (N) - L 2 (R + ,da:) : » ^ nXn 

71=1 

is an isometry. The adjoint mapping reads 

J* : L 2 (R + ,dx) - / 2 (N) : V ^ {<X«, V>»£=1- 

Set 

L(o-) = JA{p)J*. 
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L(a) is an integral operator with the kernel 

oo oo 
m=l n=l 

This can be rewritten as 

C *( x >y)=\ ifO<[y]<N 







otherwise 



Here [x] denotes the integer part of x. Notice that J* J is the identity on 
/ 2 (N) and so L(a)J = JA(a). Consequently, 

\\A(a)\\ = \\JA(a)\\ = \\L(a)J\\<\\L(a)\\. (17) 

(Hi) Denote by P n , n G Z+, the orthogonal projector onto C-Xn in -^ 2 (H^+> 

dx). 



Set 



K oS (a) = K{a) - P K(a) - K{a)P + P K(a)P - ^ P n K(a)P n . 



n=l 



In other words, we subtract from K{a) the diagonal as well as the first 
row and the first column (i.e., with index 0) with respect to the orthogonal 
system {Xn}^Lo- We can say also that the integral kernel /C° ff (x, y) vanishes 
if [x] = [y] or [x] = or [y] = and otherwise it coincides with K, a (x,y). 
Since 



P K(a)P Q -J2 P nK(a)P n 



n=l 



sup \\P n K(a)P n \\ < \\K(a) 



we have 

||X°»||< A\\K{a)\\. (18) 

(iv) It remains to estimate the norm of the difference L(a) — K oS (a). 
This is a Hermitian integral operator whose kernel does not vanish only if 
< [x] < [y] or < [y] < [x]. Suppose for definiteness that < [x] < [y\. 
Then the kernel equals, up to the multiplier —i, 



[y\ a \[y}J y 



1 \ \y 



x 



W y a J [y]-N 

1 f[y] a - [x] a y a -x a 



[y\ - [A 



y-x 



11 



Let us show that 



o <^/.(M) _!/,(£) < ^ (19) 
m \\y\J y \vj mim-m) 



First notice that 



i i f [y] a i , °{y-\y\) 

< — = -a / z-^dz < vy 11 



[vY y a J y ~ [y] a+1 

and so 

o<(— -i w-^ < * (20) 

Further set temporarily 

' " [xf y a - x a 



D 



[y] - N y - x 

a C (([*](! -t) + [y]tr l - (x(l -t)+ ytf- 1 ) dt 



The integrand in the last integral equals 

a(l-a^r 2 ((x-mi-t) + (y-[y})t) 

where £ f is a real number lying between [x](l — t) + [y]t and x{l — t) + yt. 
Notice that 

0<{x-[x)){l-t) + {y-[y))t<l. 
We assume that < o < 1. Therefore 

W' 1 - [x]*- 1 



< D < a(l - a) [ ([x](l - t) + [y]t) a ~ 2 dt = -a ■ 
Jo 



and so 

< —D < /r V . < FW1 rT -. 21 

Inequalities $2$ and ((2T|) jointly imply (fl"9l) . 

(wj From estimate (fT9|) one can deduce that L(a) — K oS (a) is a Hilbert- 
Schmidt operator and 



\\L(a)-K° s (a)\\ m <^—a. (22) 



V2 



7T 



2 
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Actually, 

POO POO 

2 dx dy\£ ff (x,y)-K?{x,y)f 

Jl J\x]+1 



poo roo 

\L(a)-K° s (a)\\ 2 s 



"OO I POO I 

2 ' - 1 ' - 1 



8ff 2 



OO 



V2 



71 



2 



vfc = l 

fwj Inequalities fTTj) . fTRjl . ifTfij) and imply that 

11^)11 < 11^)11 < H^MII + II^-^V)!! <47rV+^a. 

This shows the lemma. □ 

Lemma 3. Let A(o~) be an operator in l 2 (N) with the matrix entries in the 
standard basis 



A{a) 



for m = n 

— minjf— Y , (— Y\ for m ^ n 



Then A{a) is bounded for all < a and its norm satisfies the estimate 

\\A(*)\\<«+h£ + A\ vrV. 
Proof. Let us first show that 

|K0)||<7T. 

For cr = we get 

i 

A(0) mn = if m ^ n. 

n — m 

Considering the natural embedding Z 2 (N) C Z 2 (Z) let us denote by P + the 
orthogonal projector onto l 2 (N) in Z 2 (Z). Let B be an operator in l 2 (Z) with 
the matrix 



B mn = q{n — m ) where q{n) 



for n = 
- for n ^ 
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One can identify A(0) with P + BP + . B is a convolution operator and therefore 
it is diagonalizable by the Fourier transform T : Z 2 (Z) — > L 2 ([0, 2tt], dQ). In 
more detail, 

(TB^-^){9) = q{6)^{6) where q(0) = ^ q{n) e md . 



One finds that q{6) = —n + 6. Consequently, 

|U(0)|| = \\P+BP + \\ < \\B\\ = WFBF^W = max \q(0)\ = vr. 

6»e[0,27r] 

Suppose now that < m < n. Notice that 

w<r + 1) _xw)„ = -i Q° 

and 

(AW-xo))„=-i/.(^). 

Using Lemma [T]and Lemma [2] one can estimate 

PMII < |K0)|| + ||A(a-M)-A(0)|| + \\A(a-[a]+l)-A(a-[a] 
+ ...+ 

< 7T+ ^ + 4]7t 2 ( ( t-[ ( t])+24[ ( t] 

< 71 + | ^ + 4 ) 7T 2 CT. 

This proves the lemma. □ 

3 Boundedness of the operator i J2T=o Pk{s)Pk(s) 

We consider i Y^k=o Pk( s )Pk( s ) m the time independent frame, i.e. the oper- 
ator Q(s) defined by 

oo 

Q(s) = iV(s)* Pk{s)P k (s)V(s) = -iV(s)*V(s) = iV(s)*V(s). (23) 

k=0 

The operator V(s) is defined in (JB|). Q(s) is symmetric and its matrix entries 
in the basis {<^ n (0)} are 

(^m(O), Q(s)(f n (0)) = i (<f m {s), <fn(s)) ■ 
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Since (p n (s) depends on s only through the absolute value it holds true that 
Q(—s) = —Q(s) for s 7^ 0. For s = the operator-valued function Q(s) has 
a discontinuity. The goal of this section is to show that the operator Q(s) is 
in fact bounded. 

To compute the matrix entries one can use the identity 

ip m (s),H(s)tp n (s) 

(ipM, <p n (s)) = \ ( \ WT^ • ^ 24 ) 

Let us emphasize once more that the scalar product on the RHS should 
be interpreted as a quadratic form since, in general, <p n (s) £ Domif(s). 
The derivation goes through basically as usual even though one cannot use 
the scalar product directly. Differentiating the equation on eigenvalues one 
arrives at the equality 

H(s)(p n (s; r) + H(s)ip n (s; r) = X n (s)(p n (s; r) + A n (s)y? n (s; r), 

valid for any r > 0, in which one should substitute for H{s) and H(s) the 
corresponding formal differential operators. Next one multiplies the equality 
by np m (s; r) and integrates the both sides from e to infinity for some e > 0. 
In the integral 

poo 

(p m (s;r)d r rd r cp n (s;r) dr 



occurring on the LHS side one integrates twice by parts. Checking the asymp- 
totic behavior of the eigenfunctions near the origin, 

/ d\(M+1)/2 / 0n \ \l/2 

^)~(f) ( r( n + M + i) ) rW ( 1 + °M) 

(25) 

one finds that 

lim r(p m (s;r)d r (p n (s;r) = lim r (d r <p m (s; r)) (p n (s; r) =0. 

r— >0+ r— >0+ 

Hence sending e to actually leads to equality (j24*l) . 

Lemma 4. The matrix entries of the operator Q(s) for s ^ are given by 
the formulae 

((p m {0),Q{s)(p n (0)) = form = n, 



and 



(^m(O), Q{s)<p n [0)) = — — — r mm <^ — — , — — \ form^ n, 
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where 

7n(s) = ~ { ) ■ (26) 

Proof. Assume that m < n and s > 0. Using the explicit expression for the 
generalized Laguerre polynomials, 



n + a )L x k 
n — k ) k\ 



one finds that 

(p m (s),H(s)tp n {s)) = 2sc m (s)c n (s) 



x 



s c m( s ) c n( s ) [ ~q ) 'S'mjn 



where 



o ^ \~~V_i \ k+£ T(m + s + l)T(n + s + l)T(k + e + s) fm\ fn\ 

" T(k + s + l)T(£ + s + l)m\n\ WW' 

In this expression only the summand with k — does not vanish since 



F = for j = 0,1,..., n- 1, 



Hence 



T(m + s + l)T(n + s + l) v^,,^ T(£ + s) (n 
m ' n ~ T(s + l)m\n\ j-^ j r(£ + s + l) 

r(m + s + l)r(n + s + l) n . 

= — 77? \ . ■ -S(s,n + 1) 

r(s + l)m!n! V 7 

T(m + s + 1) 

sm! 

Furthermore, A„(s) — A m (s) = 2_B(n — m) and so 

/ ^ ^\ Y 2 Y c m {s)c n {s) T(m + s + 1) 



B ) 2B{n — m) m\ 

Now it suffices to plug in the explicit expressions for the normalization con- 
stants c m (s) and c n (s). □ 
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Using the Stirling formula one can check the asymptotic behavior of the 
matrix entries of the operator Q(s) for m and n large. It turns out that 
the operator Q(s) is in some sense close to a Hermitian operator A(s) in 
L 2 (R + ,rdr) with the matrix entries 

(<p m (0),A(s)<p n (0)) = for m = n, (27) 

and 

<*.«», A<.K(0)> = minjf^)"" 2 , fli±i N " l/2 ' 

2{n — m) \ \ n + 1 J \m + 1 

for m ^ n. (28) 

Note that -4(0+) = Q(0+). We shall also write Q(s) mn instead of 
(¥>m(0), Q(s)<p n (0)), and similarly for A(s). 

Lemma 5. Let A(s) be the Hermitian operator in L 2 (EL +J rdr) defined by 
relations f21\) and (H2Jj. Then Q(s) —A(s) is a Hilbert- Schmidt operator and 
it holds true that 

||Q( s )-A( s )||hs<^|(1 + |,|)( 3+ '^ 2 . 

Proof. Let us suppose for definiteness that s > and m < n. For x > 1 set 

T(x + s) 



9s(x) 



x s T(x) 



One can express 



\Q(s) mn - A(s) mn \ = — — — \g s (m + l) 1 ' 2 - g s (n + l) l ' 2 \ 



m + 1 x 



71 + 1 



^(n + 1)- 1 / 2 



< 



1 Z" 1 

- g s (n + I)" 1 / 2 / g s (m + 1 + (n - m)t)^ 1/2 



Notice that 



4 

x |^g(m + 1 + (n — m)t) \ dt. 
g' s (x) r'(x + s) V(x) s 



g s {x) F(x + s) T(x) x 

Using the well known formula for the logarithmic derivative of the gamma 
function, 

r'(*) 

T(z) 



I + 7 + f;(_5__I), 

z ^— ' V w + z n J 

n=l ' 



(29) 
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one finds that 



g s (x) \ ^ (n + x){n + x + s) x 



v n=0 

oo 

< S ' ' 



< S 



™ 1 f 

^ (n + x) 2 x 

C 2 Jx 



n=0 

x" j x y 2 x 



x 



Similarly, 



9s(x) " \J X y(y + s) x 

s 

X 



> 



s 2 



2x 2 



In particular, 

W s (x)\< S -^g s (x). 

Jb 

From here one derives the estimates, for t G [0, 1], 

g s (m+l + (n-m)t) _ /_ T+ 1 gjjy) 



V >/ m+l+(n—m)t 9s 

(v) 

Wi \y \ y 

= exp ( (m + 1 + s) In ( 1 + 



< cxp( / (- - \n(^—^ ] I d/y 



m + 1 



(n + l + s)ln(l + 



V n + 1 



< (l + s) 1+s 



and 



4 + 1)V2 y ( m + ! + („ _ m ) t )2 

df 



^ r (1+s)(3+s)/2 fr^7- — 

4 y (m + 1 



+ (n — m)t)" 

18 



Let F(t) be a Hermitian operator in L 2 (M + ,rdr) with the following matrix 
entries in the basis {y? n (0)}: 

F{t) mn = for m = n, 

and 

F(t) mn = (m + 1 + (n — m)t)~ 2 for m < n. 
Then F(t) is a Hilbert-Schmidt operator and 

oo oo 

= 2 E E (m + l + (n-m)t)- 4 



m=0 n=m+l 

oo 



m=0 



(m + 1 + ft/) 4 
E (m+ l) 3 



< 
~ t 



2 

m=0 

1 



Hence 



s)-A( S )|| HS < ^(l + «) (3+s)/2 ^ IW)llHsdt 

< 1 (l + s )(3+.)/2_ 

This proves the lemma. □ 

Combining Lemma El and Lemma El we deduce that the operator Q(s) is 
actually bounded. 

Lemma 6. The operator Q(s) is bounded and its norm satisfies the estimate 

\\Q(s)\\<l + 12\s\ + l\s\(l + \s\)^ 2 . 

Proof. Let A(s) be the Hermitian operator in L 2 (M + ,rdr) defined by rela- 
tions (J2TJ) and (|28j). According to Lemma El it holds true that 

IW'>IIS5(" + (t + 4 Vt 
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Lemma El leads to the estimate 

\\Q(s)\\ < \\A(s)\\ + \\Q(s) - A(s)\\ 

< ^(-+(^ f + 2) 7 r 2 | S | + | S |(l + | S |)^ 2 ). 

Since (l + 1/(6^)) vr 2 < 12 the lemma follows. □ 

4 The meaning of the propagator U T (s, sq) 

As already discussed in the Introduction the natural propagator U t (s,sq) 
defined in (JHJ is not related in the standard way to the Hamiltonian tH(s) 
defined in (QJ. In particular it is not clear if U t (s,Sq) maps the domain 
Domif(so) into Domif(s). This is why we propose in the Appendix the no- 
tion of a propagator weakly associated to a Hamiltonian, see Definition IA.31 
We should like to emphasize that this relationship is unique, i.e. at most one 
propagator can be weakly associated to a Hamiltonian. 

In this section we show that U T is weakly associated to tH and that 
(s,r) i— > U T (s, s )ijjo(r) satisfies the Schrodinger equation as a distribution 
for all Vo e L 2 (IR + , rdrdip). 

Proposition 7. The propagator U T (s, Sq) is weakly associated to tH(s). 
Proof. Relation (jSJ means that 

U T (s,s ) = V(s)e-^C(s,s )e^ n ^V(s )- 1 . 

So starting from C(s, So) one can reach U T (s, Sq) by two consecutive unitary 
transformations. The propagator C(s,s ) was defined in Q. It corresponds 
to the Hamiltonian — Q T (s) defined in 0. According to Lemma El the func- 
tion ||Qt(s)|| = ||Q( S )|| is locally bounded and thus C(s,sq) is given by the 
Dyson formula, see relation (p?T]) in Section El 

First we apply Proposition IA.4I in which we set A(t) = —Q T (t), D = 
Dom#(0), T(t) = exp(-irfi(t)) and 

X(t) = % (d t e- irm ) e lTm = rW(t). 

We conclude that the propagator e~ tTn( - s ' > C(s, so)e"~ n(so ) is weakly associated 
to 

tW{s) - e- iTn{s) Q T (s)e iTn{s) = tW{s) - Q(s). 
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Next we apply Proposition IA.6I in which we set H(t) = rW(t) — Q(t) and 
U(t,s) = e- ir ^C{t,s)e ir ^ s \ Recall further that V(t) was defined in ©. 
We conclude that U T (s, sq) = V(s)U(s, so)V(so)" 1 is weakly associated to 

tV(s)W(s)V(s)- 1 - V(s)Q(s)V(s)- 1 + z\>(s)\/(s)" 1 = tH(s). 

The proposition is proven. □ 
In the studied model Tt = L 2 (M + , rdr) and so 

K = L 2 (R, H, ds) = L 2 (Rx R+, rdsdr). 

Let S) = H(s) ds be the direct integral of the family of self-adjoint opera- 
tors H(s) which is nothing but a multiplication operator in /C. Let K T be the 
quasi-energy operator associated to the propagator U T (s, s ) (see Appendix). 
According to Proposition it holds true that 



K T = -id s + tS) . (30) 

To an initial condition ipo G Ji we relate the function 
i/j(s,r) = (U T (s,0)ip )(r) which is a locally square integrable function in the 
variables s and r. We now show that ip(s, r) fulfills the Schrodinger equation 
in the space of distributions f (Ix ]0, oof). Let us note that for the proof 
it suffices to know that — id s + tS) C K t , the stronger property (TTTTj) is not 
necessary. 

Proposition 8. For every tpo £ H, the function ^(s^r) = (U T (s, 0)i/; )(r) 
satisfies the Schrodinger equation in the sense of distributions. 

Proof. Let £ G Cq°(Mx ]0, +oo[) be an arbitrary real-valued test function. 
Set g(s,r) = £(s,r)/r. Clearly, g G Dom(— id s + rfj) C DomK T . Let 
[a, b] x [c,d] be a rectangle containing supp£ and choose r] G C£°(R) so that 
T) = 1 on a neighborhood of the interval [a, b]. From Proposition lA.2l we know 
that K T (rj(s)ip(s ) r)) = —ir]'(s)ip(s,r). From the choice of rj it follows that 

= -i(g,r)'4>))C = {g,K T {jYil))) K = ((-id s + rfyg^t/j)^ 

The last term equals 

id s - r) + tH(s)- £(s, r) Jt7(s)'0(s, r) rdsdr 
r r I 



id s £(s, r) + r ( — <9 r r<9 r — I — - ( s H ) I £(s, r) ) t/>(s, r) dsdr. 
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This means that 



1 „ „ If Br 2 ^ 



-id s ip(s, r) + t\ -- d r rd r + — + -y-J jip(s, r) = 
in the domain Ix ]0, +oo[ in the sense of distributions. □ 



5 Proof of the adiabatic theorem 

We follow the strategy explained in the Introduction. The adiabatic propa- 
gator Uad (see (jlj)) and the propagator U T defined in (jSJ differ by C defined 
by ®. Since Q T (s) = e iT ^ {s) Q(s)e- irn{s \ defined in 0, is unitarily equiva- 
lent to Q(s) it is bounded, uniformly in s on every bounded interval [0,5]. 
Hence C(s, s ) exists and is given by the Dyson formula: 

C(s,s Q ) = l + ^2i n / dst ds 2 ... ds n g r (s 1 )g r (s 2 ) . . .Q T (s n ). 

n=l ^ s o J sa Jsq 

(31) 

The task is to estimate the norm of the integral of Q T . This will be done 
by the integration by parts technique developed in the following two lemmas. 

The first step is to find a bounded differentiable solution X(s) of the 
commutation equation 

Q(s)=i[W(s),X(s)}. 

The operator W(s) was defined in Q. The off-diagonal entries of the X(s) 
are determined unambiguously, 

/ fn\ vf \ (n\\ ■ (^m(0),Q(g)^»(0)) . . 

(cp m {0),X(s)ip n (Q)) = -i — —— — - — (32) 

A m [S) — A n [S) 
Sgn(s) . (j m (s) 7 n (s) \ 

mm < — — , — — > for m fi n, 



AB(n - m) 2 [ 7„(s) ' 7 m (s) 

with 7n(s) defined in (|26| . We set 

(<p m (0),X(s)(p n (0)) = for m = n, (33) 

and write again X(s) mn instead of (ip m (0),X(s)ip n (0)). 

Lemma 9. The operator X(s) defined by relations pity) and pJOfy is bounded 
and its norm satisfies the estimate 



\X(s)\\ < 



TV 2 



12 B 
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The derivative X(s) exists in the operator norm and satisfies the estimate 

(1 + V2)n 2 



\\X(s)\\ < 



48 B 



Proof. The operator norm of X(s) is bounded from above by the Shur- 
Holmgren norm, 



1 00 1 

1 x — - 1 7T 



||X( S )||<||X( S )||sh= sup J2\X(s) mn \<— 



2 



2B ^ k 2 12 B 

k=l 



Suppose that s > and m < n. Let us estimate the derivative of X(s) 
Using one finds that 

7m(s)Y 7m(s) fT'(m + s + l) T'{n + s + l) 



ln(s)J 2 7n (s) Vr(m + s + l) T(n + s + l) 

- m (s) ^ " - ni 



2 7n(s) (k + m + s + l)(fc + n + s + 1) 



Hence 

dy 



< ^- - ______ + 



oo 



8B(n — m) \(m+l)(n + l) J 1 (y + m){y + n) 

1 . f n+l 
H In 



8B(n — m) \(m+l)(n + l) n — m \m+l 
Thus we get, for m / n, 



ds 



<J_f 1 + 1 ). (34) 

8B \{m + l)(n + 1) |n — m| min{m + 1, n + 1} 



Let be a Hermitian operator in L 2 (M + ,rdr) with the matrix entries 

dX(s) mn /ds. From the estimate (|34"|) we deduce that X(s) is a Hilbert- 
Schmidt operator and 



/ oo oo 



1/2 



85 \ ^ (m + l) 2 ^ (n + l) 2 

\m=0 v 7 n=0 v y 



oo ^ oo 



1 ' " I " 1 1/2 

+ 8B [ 2 S (m+ l) 2 S i (n_ m ) 



m=0 
(l + V^TT 2 

48 5 
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Furthermore, since estimate (|34*|) is uniform in s one can apply the Lebesgue 
dominated convergence theorem to conclude that 



lim 



1 



(X(s + e)-X(s))-X(s) 



0. 



HS 



Hence the derivative of the operator-valued function X(s) exists in the op- 
erator norm and equals X(s). □ 

The matrix entries of the operator Q T (s) defined in (0) equal 

(<p m (0), Q T (s)ip n {0)) = % e^W^W (ip m (s), <p n (s)) . 

Notice that the both operators Q(s) and W(s) = Q'(s) are diagonal in the 
basis {(^n(0)} and therefore they commute. 



Lemma 10. It holds true that 

Q T (u) du 



<n + l±^i.i 



7T 



6Br 



Proof. Suppose that s > 0. The integral can be rewritten as follows, 

Q T (u)du = i f e iTn{u) [W(u),X(u)}e- iTn{u) du 
o Jo 

= I J' ((e^M)' X(u) e- trn ^ + e iTn{u) X(u) ( e -^(«))') du 

= \ J ({e lrn(u) X(u) e - irn{u) )' - e lrn{u) X(u)e- iTn{u) ^j du . 

Consequently, 

jT Q T (u) du = i (V Q(s) X(s) e-^ (s) - X(0) 

-J e iTniu) X(u)e- iTn(u) d u y 

More precisely, the derivation of this equality was rather formal but it 
becomes rigorous when sandwiching the both sides with the scalar product 
(^m(O), ■ (p n (0)). This is to say that the both sides have the same matrix en- 
tries in the basis {(p n (0)}. But since the equality concerns bounded operators 
it holds true. 
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Using Lemma El one arrives at the estimate 

1 



Q T (u) du 



< 



< 



71 



*( a )|| + pr(o)|| 

1 1 + V2 



\X(u)\\du 



Bt \ 6 



48 



The lemma is proven. 



□ 



We can now show that the adiabatic propagator Uad(s,0) (see (J2J)) is 
close to the propagator U T (s, 0) = Uad(s, 0)C(s, 0) defined in provided 
the adiabatic parameter r is large. 

Proposition 11. It holds true that 



\U T (s,0)-U AD (s,0)\\<M(s)e^ M ^ 



71 



3Bt 



where 



M(s) = ^ + 12\s\ + ^\s\(l + \s\f + ^ 2 . 



(35) 



Proof. According to Lemma El < M(s), and from Lemma HH one 

easily deduces that 



Q T {u) du 

Using formula JSJ) one can estimate 
||C/ r ( S ,0)-^( S ,0)|| = ||C(s,0)-I|| 

°° f\s\ 



< 



7T 



3fir 



M{s). 



< 



n=l 

X 

7T 



ds n _i ||Q T (si)|| . . . ||Q T (s n _i 



o 



71 

3Bt 



Ym( 

3Bt ^ v 

n=l 

oo 

n=l 

The proposition is proven. 



ds n Qr^nj 



*)"' I dsi . 
o 

(n-l)! ' 



ds ra _! 



□ 



25 



6 The general dependence on time 

Here we show that the adiabatic theorem extends to Hamiltonians of the 
form 

H<(s) = H(((s)) 

where H(s) is defined in (0) and ( G C 2 (IR) is a real- valued function. In order 
to simplify the discussion and to avoid considering discontinuities (recall that 
Q(s) is discontinuous at s — 0) we shall further assume that C'(s) > and 

C(o) = o. 

Set 

V c (s) = V(((s)), W<(s) = W(((s)), &(s)= f S W c (u) du. 

Jo 

Let C^(s, so) be the propagator related via the Dyson formula to the Hamil- 
tonian — Qf(s) where 

Qi(s) = exp(zrO^( S )) Q<(s) exp(-zrO<( S )) , Q<(s) = C(s)Q(C(s)). 

Exactly in the same way as in the proof of Proposition one can show that 
the propagator 

U<(s,s ) = V^s) exp(-zYO<^)) C^s,s )exp(ir^(s Q )) V^so)' 1 

is weakly associated to the Hamiltonian H^(s). The adiabatic propagator 
now reads 

U c AD (s, s ) = V<(s) exp(-ir(fi c (a) - &(s ))) V<(s )-\ 

Proposition 12. Assume that ( G C 2 (R), ('(s) > and C(0) = 0. Then 
there exists a locally bounded function m^(s) such that 

m^(s) 



Vs G 



U^s,0)-Ui D (s,0) 



< 



Bt 



Proof. Suppose for definiteness that s > 0. Recall that < M(s) 

where M(s) was defined in The operator- valued function 

X C (s) = C(s)X(((s)), 
with X(s) being defined in (E?2*|l and ijHHJl . satisfies the commutation equation 

Q^s)=t[W^s),X<(s)}. 
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Quite analogously as in the proof of Lemma [TQl one derives the estimate 



Q C T (u) du 



In virtue of Lemma El we have 



<- [\\X<(3)\\ + ||^(0)||+ / \\X<(u)\\du . 



\\xH*)\\<£bC{*) 



and 



\\X c (u)\\du < 



7T 



2 rs 



12B 



\C'(u)\du + 



:i + v / 2) 

485 



7T 



2 /-s 



Civ) 'du. 



Hence 



where 



Qi(u) du 



< 



Br 



<7 C (*) = ^ I C(0) + sup C'(«) + / |C»I dti + 



12 



0<-u<s 



1 + ^2 



C(uydu 



Finally one can proceed similarly as in the proof of Proposition ITT1 to derive 
the estimate 



( r c{s) 

U<(s,0)-U c AD (s,0) = \\C<(s,0)-l\\ <expf J 



M(v) dv 



Br 



This completes the proof. 



□ 



Appendix. Propagator weakly associated to a 
Hamiltonian 

By a propagator U(t, s) we mean a family of unitary operators in a separable 
Hilbert space TC depending on t, s e K which satisfies the conditions: 

(i) U(t, s) is strongly continuous jointly in t, s, 

(ii) the Chapman-Kolmogorov equality is satisfied, i.e. 

Vt, s,r el, U(t, r)U{r, s) = U(t, s). 
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Let H(t), t G R, be a family of self-adjoint operators in H. The domain may 
depend on t. The standard way how one relates a propagator U(t, s) to H(t) 
is based on the following two requirements: 

(i) Vt,seR, U{t,s)(DomH(s)) = Dom H{t), 

(ii) G Dom#(s),Vi G R, idtU(t,s)ip = H{t)U{t, s)ip. 

Clearly, if a propagator exists then it is unique. In some situations, however, 
these requirements may turn out to be unnecessarily strong. In particular 
this is true for the model studied in the current paper. The heart of the 
problem is illustrated on the following example. 

Let A(t) be a family of bounded Hermitian operators in 7i which is uni- 
formly bounded. Then the propagator exits and is given by the Dyson 
formula. Let us call it C(t,s). Let D C TC be a dense linear subspace, 
and let T{t) be a strongly continuous family of unitary operators such that 
D is invariant with respect to T{t) and for every ip G D there exists the 
derivative d t T(t)ip. Furthermore, suppose that X{t) = iT(t)T(t) _1 , with 
DomJ(t) = D, is a self-adjoint operator for all t (the dot designates the 
derivative). A formal computation gives 

T(t) (-id t + Ait^Tit)' 1 = -id t + X(t) + T{t)A{t)T(t)- 1 . 

If C(t,s) preserved the domain D then the propagator T(t)C(t, s)T(s) -1 
would solve the Schrodinger equation for X(t) +T(t)A(t)T(t)~ 1 on D. Thus 
it is natural to associate it to this family of self-adjoint operators. The 
hypothesis on C(t, s) need not be, however, satisfied since A(t) is an arbitrary 
family of bounded operators and so C(t,s) will in general not preserve this 
domain. 

In this appendix we propose a way how to associate a propagator to 
a given time-dependent Hamiltonian in a weak sense. This association is 
more general than the standard one (which supposes a constant domain and 
solving the Schrodinger equation in the strong sense) and it is still unique 
(i.e.: there is at most one propagator weakly associated to a given time 
dependent Hamiltonian). 

Here we develop this approach only to an extent which makes it possible 
to apply these ideas to the studied model with a time-dependent Aharonov- 
Bohm flux. In particular, the described example is covered by Proposi- 
tion EU below. 

Let X be a Banach space. We shall say that a vector-valued function 
/ : R — > X is absolutely continuous on R if it is absolutely continuous on 
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every compact interval / C R. By the symbol AC (til, X) (or just AC if 
there is no danger of misunderstanding) we shall denote the space of all 
absolutely continuous vector- valued functions f(t) such that the derivative 
/'(£) exists almost everywhere on til. In such the function is 

locally integrable and f(t) = f(0) + /*/'(«) ds HU Theorem 3.8.6]. If the 
Banach space A? has the Radon-Nikodym property then the space AC (til, X) 
coincides with the space of absolutely continuous vector-valued functions 
AC (til, X). Let us recall that X is said to have the Radon-Nikodym property 
if the fundamental theorem of calculus holds, i.e. if any absolutely continuous 
function is the antiderivative of a Bochner integrable function. For example, 
separable Hilbert spaces are known to have the Radon-Nikodym property [|]. 

Clearly, if /, g G AC(til, TC) then the function (f(t),g(t)) is absolutely 
continuous and 

dt(f(t),g(t)) = {fit), g{t)) + (f(t),g'(t)) a.e. 

Similarly, if A G AC (til, 13(H)) and / G AC(R,H) then A(t)f(t) G AC(R,H) 
and 

d t A(t)f(t)=A(t)f(t) + A(t)f'(t) a.e. 

Let {efc} be an orthonormal basis in Ti. A vector-valued function f(t) = 
J2Vk(t)^k belongs to AC(til,TC) if and only if the following two conditions 
are satisfied: 

(i) 3a G til such that \rjk(a)\ 2 < oo, 

(ii) V*, Vk G AC, and (£, W k (t)\ 2 f" G L\ ac (til) . 

From here one easily derives the following criterion (alternatively, one can 
again consult [101 Theorem 3.8.6]). 

Lemma A.l. A vector-valued function f : til — > Tt belongs to AC (til, H) if 
and only if the following two conditions are satisfied: 

(i) there exists a total set T C 7~t such that for all ip G T, (ip,f(t)) is 
absolutely continuous, 

(ii) the derivative f'(t) exists a.e. and G Li oc (R). 

Set K = L 2 (til, Ti, dt). Let us recall that to every propagator U(t, s) on Ti 
one can relate a unique self-adjoint operator K in /C which is the generator 
of the one-parameter group of unitary operators exp(— iaK), a G R, defined 
by 

(e-™ K f)(t) = U(t,t-a)f(t-a) 
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[TT| . K is called the quasi-energy operator. Equivalently, 

/■© 

K = il(-id t )if where U= / C/(t, 0) dt. (A.l) 



So / G DomX if and only if U(t, 0) G Dom(— z"9 t ) which means that 

/ G L 2 , C/(t,0)-7(0 G AC and (C/(t, O)- 1 /^))' G L 2 . 

From (|A.lj) one concludes that the spectrum of .K" is purely absolutely 
continuous and coincides with R. So the kernel of K is always trivial. It 
seems to be natural, however, to introduce a generalized kernel of K, called 
Ker K , as follows: 

KeroiT = {/ G L 2 0C (M, ft, dt); V77 G C °°(R), 77/ G Domif 
and K(r)f) = -irj'f}- 

Since K can be very roughly imagined as the formal operator —id t + H(t) the 
elements of Ker K can be regarded as solutions of the Schrodinger equation 
in a weak sense. 

Proposition A. 2. Let U(t, s) be a propagator and let K be the quasi-energy 
operator associated to it. Then it holds 

Ker K = {C/(t,0)V; tpeH}. 

Proof. If f(t) = U{t,0)ip, with ip G H, and rj G C£°(R) then, in /C, there 
exists the derivative 



d 



d 



i ^ {e- laK vf)(t) ^ = 1— {v(t - a)U(t, Oty) n = -irf{t)f{t). 



(7=0 dcr 



<7=0 



Hence, by the Stone theorem, 77/ G Dom.fr and K(rjf) = —irj'f. 

Conversely, suppose that / G Ker A" and set g{t) = U(t, 0) _1 /(t). Let 77 
be a test function. From 1JA.1J) one deduces that 77(/ G Dom(— and 



dt{vit)g{t))= V '(t)g(t) 



a.e. 



Since 77 G C^°(R) is arbitrary this implies that g(t) G AC(R, H) and g'(t) = 
a.e. Consequently, g(t) = ip G 7i is a constant vector-valued function and 

/(t) = [/(t,o)v>/ □ 

It is known that the correspondence between the propagators and the 
quasi-energy operators is one-to-one [Til Remark (1) on p. 321]. On one hand, 
by the very definition, K is unambiguously determined by U(t,s). On the 
other hand, if U(t, s) and U\(t, s) are two propagators with equal quasi-energy 
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operators, K = Ki, then U(t,s) = Ui(t,s). This uniqueness result is also 
a straightforward corollary of Proposition IA.21 Actually, Proposition IA.2I 
implies that for every if) G Ti there exists ipi G Ti such that U(t, 0)if> = 
U\(t,0)ipi for all t (we use the strong continuity of the propagators). By 
setting t = one finds that if) = if)\. Hence U(t,0)if> = U\(t,0)ip for all 
if) ETi. Consequently, 

U(t,s) = U(t,0)U(s, 0) _1 = U x (t, 0)C/i(s, O)" 1 = C/i(t,s). 

For a family of self-adjoint operators H(t), t G R, set Sj = H(t)dt. 
This means that / G /C belongs to Domi^ if and only if /(£) G Domif(t) a.e. 
and \\H(t)f(t) \\ G L 2 {R, dt). Then ^ is a self-adjoint operator in K. In what 
follows we shall always suppose that the intersection Dom(— id t ) H Dom^ is 
dense in K. For example, this is true in the case when the domain Domif (t) 
is independent of t. Consequently, — id t + f) is a densely defined symmetric 
operator. 

Definition A. 3. We shall say that a propagator U(t, s) is weakly associated 

to H(t) if 

K = -td t + £ • (A.2) 

Notice that equality (|A.2|) is equivalent to the following two conditions: 

(i) -id t + ^CK, 

(ii) — id t + is essentially self-adjoint. 

Furthermore, it is important to note that this definition still guarantees the 
uniqueness, i.e. to H(t) one can weakly associate at most one propaga- 
tor U(t,s). Actually, if U(t,s) and Ui(t, s) are weakly associated to H{t) 
then K = Ki according to equality (IA.2J) . But due to the one-to-one corre- 
spondence between the propagators and the quasi-energy operators we have 
U(t,s) = U 1 (t,s). 

Now we are ready to formulate and prove two propositions which are 
directly applicable to the model studied in this paper. 

Proposition A. 4. Let A(t) be a family of bounded self-adjoint operators 
in 7i which is locally bounded. Let C(t, s) be the propagator associated to 
A(t) via the Dyson formula. Let D C Ti be a dense linear subspace and let 
T{t) be a strongly continuous family of unitary operators in Ti obeying the 
conditions: 

(i)Vte E, T(t)D = D, 
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(ii) G D, T(t)ip is continuously differentiable, 

(lit) Vi G R, X(t) = it{t)T{t)- 1 , with DomX(t) = D , is a self-adjoint 
operator. 

Then the propagator T(t)C(t, s)T(s) -1 is weakly associated to the family 
X{t)+T{t)A{t)T{t)-\ 

Proof. Set 

Y(t) = X(t) +T(t)A(t)T(t)- 1 , 2) =/ Y{t)dt,% = T(t)dt. 

Jr Jr 

Let K Y be the quasi-energy operator associated to the propagator T(t)C(t, s)T(s)~ 1 . 
Set 

C(t) = C(t,0), £ = / C{t)dt. 

Jr 

C(t) is a family of unitary operators which satisfies C{t) G AC(R, B(TC)) and 
A(t) = iC r (t)C(t)- 1 . 
(i) Let us verify that 

-id t + Z)cK Y = ZC(-id t )€T 1 T' 1 . 

Suppose that a vector-valued function / : R — > belongs to Dom(— z<9t + 2)). 
This happens if and only if / obeys the conditions: / G L 2 , / G AC, /' G L 2 , 
/(t) G D a.e. and Y(t)f(t) G L 2 . In that case the function r(t) -1 /(t) is 
differentiable a.e. and the derivative 

(T(t)-if(t))' = T(t)- 1 (f'(t)+iX(t)f(t)) 

is square integrable. Moreover, if tp G D then the function (ip , T (ty 1 f (t)) = 
{T(t)ip, f{t)) is absolutely continuous. According to Lemma lA. II this implies 
that T(t)-7(t) G AC(R,W) and consequently C^T^)' 1 f(t) G AC as 
well. Furthermore, a straightforward computation yields 



Y(t)f(t) = i^Tity'f^+T^C^C^Tity^it) 

= limc^ycity'nty'fit) 

= if\t)-iT{t)c{t)(c{ty l T{ty l f{t))' . 

Hence (C^Tity 1 f [t))' G L 2 , f G Domi^y and + Y(t)f(t) 

Kyf(t). 
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(ii) Let us verify that — id t + 2) is essentially self-adjoint. Suppose that 
g G Dom(— i<9 t + 2))* satisfies (— id t + 2})*g = zg with Im(z) 7^ 0. This means 
that 

V/ G Dom(-zd t + 2J), + 2))/, <?>*; = q)k- 

Choose fit) = rj{t)T{t)ip where ip G D and r\ G Co°(R) is real-valued. Then 
/ G Dom(— z<9 4 + 2J) and an easy computation shows that 

(-idt + %))f(t) = -i V '{t)T(t)if; + r,{t)T(t)A(t)if;. 

Hence for all r\ G Cq°(R) we have 

(ir/{t)(T{t)ij>,g(t)) + V (t){T{t)A{t)iP,g{t)))dt = z [ r)(i)(T(t)ip, g(t))dt . 

Jr 

Setting 

F(t) = (T(t)iP,g{t)), G(t) = (T(t)A(t)iP,g(t)), 

we find that 

-id t F(t) + G{t) = zF(t) (A.3) 

in the sense of distributions. Since both F(t) and G(t) are locally integrable, 
a standard result from the theory of distributions tells us that F(t) is abso- 
lutely continuous and equality (I A. 3 J) holds true in the usual sense. Moreover, 
equality (|A.3|I implies that 

d t (e 2Im{2) *|F(t)| 2 ) =2e 2Im{z)t lm(F(t)G(t)). 

Let us now choose an orthonormal basis {ip k } whose elements all belong 
to the domain D. Let us write F k instead of F and G k instead of G when 
replacing if) by ^fc- We have derived the equality 

\F k {t)\ 2 = e~ 21m{z)it - a) \F k {a)\ 2 + 2 f\- 21miz ^- s hm(F\^)G k {s))ds (A.4) 

J a 

which is valid for all k and all a, t G R. Observe that 
El F ^)| 2 =ll^)l| 2 a.e., 

k 

< \\g(s)\\ P( S )T( S )-^( S )|| GLj oc (R,d S ) a.e., 

and 

^i^)G fc (s) = <( ? ( S ),T( S )A( S )T( S )^ 1 ^( S )) G R a.e. 
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Summing in k in equality (|A.4J) we find that 

|^(t)||=e- Im ^)||^(a)|| 

for almost all a, t 6 R. Since ||<?(t)|| is square integrable this is possible only 
if = a.e. □ 

Proposition IA.4I has a corollary justifying the adverb "weakly" in Defini- 
tion El 



Corollary A. 5. Assume that a propagator U(t,to) is associated as a strong 
solution of the Schrddinger equation to a time- dependent Hamiltonian H{t) 
which has, however, a time-independent domain (i.e. the relationship between 
the propagator and the Hamiltonian is the usual one). Then U(t, t ) is weakly 
associated to H(t). 

Proof. In Proposition EH it suffices to set D = DomiJ(O), T(t) = U(t,0) 
and A(t) = 0. Then X{t) = H(t), C{t,s) = I and T(t)C(t, s)T(s)- 1 = 
U(t,s). □ 

Proposition A. 6. Suppose that V(t), t Gl, is a family of unitary operators 
which is continuously differentiate in the strong sense. Let H(t) , t e E, be 
a family of self-adjoint operators such that Dom H{t) = D for all tel. Set 

H(t) = VtyHtyVtt)- 1 + tV(t)V(t)-\ 

If the propagator U (t, s) is weakly associated to H{t) then the propagator 
U(t,s) = V(t)U(t, s)K(s)" 1 is weakly associated to H(t). 

Proof. Set 

/■© /■© 
U(t) = U{t,0), 11= / U(t)dt,V3= V(t)dt. 



By the assumption, iX(— id t )ii 1 = — id t + fj. We have to show that 

93H(-i8| t )ir 1 ar 1 = -%d t + 53 . 

Since 



it is sufficient to verify that 

w{-id t + s))^- 1 = -id t + ft. 
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This would also imply that Dom(— id t ) PI Dom(^i) is dense in /C. 

A vector-valued function / : R — > 7i belongs to Dom(2J(— i9 t + i^)53 _1 ) if 
and only if it satisfies the conditions: / G L 2 , ^(t)- 1 /^) e AC, e 
L 2 , V{t)- l f{t) E D a.e. and H{t)V {ty 1 f{t) E L 2 . Let us note that from the 
continuous differentiability of V(t) in the strong sense and from the uniform 
boundedness principle it follows that V(t), t E M, is a family of bounded 
operators which is locally bounded. Furthermore, V(t)* = V{t)~ l is contin- 
uously differentiable in the strong sense as well and Vit)^ 1 ^ E AC for all 
ip E H. Suppose that / E L 2 . If V^(t) _1 /(*) e AC then /'(*) exists a - e - and 
is locally integrable, the function {ip, f(t)) = (^(t)- 1 ^, V(t)"7(t)) is 
absolutely continuous for all ip E TC and therefore, by Lemma lA - !} f(t) E AC. 
Similarly, the converse is also true. If f(t) E AC then V{t)~ l f{t) E AC. 

Using these facts and the relation between H(t) and H(t) (including that 
Dom H(t) = V(t)D) one easily finds that the domains of %3(—id t + fj)2J -1 
and — id t + ft coincide and that 

V(t) {-id t + Hit^Vity'fit) = -if(t) + H(t)f(t) 

for every / E Dom(—id t + $)). □ 

Remark. Proposition IA.6I can be easily extended to the case when the family 
of unitary operators V(t) is continuous and piece- wise continuously differen- 
tiable in the strong sense and in each point of discontinuity there exist the 
limits of the derivative both from the left and from the right. 
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